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Effects of Imperfect Gate Operations in Shor's Prime Factorization Algorithm 
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The effects of imperfect gate operations in implementation of Shor's prime factorization algorithm 
are investigated. The gate imperfections may be classified into three categories: the systematic error, 
the random error, and the one with combined errors. It is found that Shor's algorithm is robust 
against the systematic errors but is vulnerable to the random errors. Error threshold is given to the 
algorithm for a given number N to be factorized. 
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I. INTRODUCTION 

Shor's factorization algorithm [l| is a very impor- 
tant quantum algorithm, through which one has demon- 
strated the power of quantum computers. It has greatly 
promoted the worldwide research in quantum computing 
over the past few years. In practice, however, quantum 
systems are subject to influence of environment, and in 
addition, quantum gate operations are often imperfect 
0, Q . Environment influence on the system can cause de- 
coherence of quantum states, and gate imperfection leads 
to errors in quantum computing. Thanks to Shor's an- 
other important work, in which he showed that quantum 
error correlation can be corrected [4J. With quantum 
error correction scheme, errors arising from both deco- 
herence and imperfection can be corrected. 

There have been several works on the effects of deco- 
herence on Shor's algorithm. Sun et al. discussed the 
effect of decoherence on the algorithm by modeling the 
environment @ . Palma studied the effects of both deco- 
herence and gate imperfection in ion trap quantum com- 
puters [6j. There have also been many other studies on 
the quantum algorithm 0, H, H, 03 • 

The error correction scheme uses available resources. 
Thus it is important to study the robustness of the algo- 
rithm itself so that one can strike a balance between the 
amount of quantum error correction and the amount of 
qubits available. In this paper, we investigate the effects 
of gate imperfection on the efficiency of Shor's factoriza- 
tion algorithm. The results may guide us in practice to 
suppress deliberately those errors that influence the algo- 
rithm most sensitively. For those errors that do not affect 
the algorithm very much, we may ignore them as a good 
approximation. In addition, study of the robustness of 
algorithm to errors is important where one can not apply 
the quantum error correction at all, for instance, in cases 
that there are not enough qubits available. 



The paper is organized as follows. Section II is devoted 
to an outline of Shor's algorithm and different error's 
modes. In Section III, we present the results. Finally, a 
short summary is given in Section IV. 



II. SHOR'S ALGORITHM AND ERROR'S 
MODES 

Shor's algorithm consists of the following steps: 
1) preparing a superposition of evenly distributed 
states 
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where q = 2 L and N 2 < q < 2N 2 with N being the 
number to be factorized; 

2) implementing ?/ a modiV and putting the results into 
the 2nd register 
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3) making a measument on the 2nd register; The state 
of the register is then 
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where j < 

4) performing discrete Fourier transformation (DFT) 
on the first register |</> 3 ) = [J2 c f (c) |c) ) \z), where 
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This term is nonzero only when c = k-, with k — 
0, 1, 2...r — 1, which correspond to the peaks of the dis- 
tribution in the measured results, and thus this term be- 



comes /(c) = -^=e i . The Fourier transformation is 
important because it makes the state in the first register 
the same for all possible values in the 2nd register. The 
DFT is constructed by two basic gate operations: the 

single bit gate operation Aj = ^ j "1 ) ' w ^ c ^ * s 
also called the Walsh-Hadmard transformation, and the 
2-bits controlled rotation 
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with 6jk 
DFT is 



The gate sequence for implementing 



(A q -l)(B q -2q-lA q -2) ■ ■ ■ (Boq-lBoq-2 ■ ■ ■ BqiAq). 

Errors can occur in both Aj and Bjk- Aj is actually a 
rotation about y-axis through ^ 



Aj{0) = e* b * v = /cos(-)-»sin(-)a y 



cos(|) -sin(f) 
sin(f) cos(f) 



If the gate operation is not perfect, the rotation is not 
exactly ^. In this case, Aj is a rotation of ^ + 25 



cos(<5) — sin(<5) — (sin(5) + cos(S)) 



y/2 V sin(5) + cos((5) cos(<5) — sin(<5) 
If 6 is very small, we have: 

3{) ~^2\l + S+ 1-5 

Similarly, errors in Bjk can be written as 
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With these errors, the DFT becomes 
l«> - 4= E e^ +S ' )a (l + S' C )\S) 
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where 5 C and 5' c denote the error of Aj and Bjk, respec- 
tively. 

Let us assume the following error modes: 1) system- 
atic errors, where 5 C or 5' c in (1) can only have system- 
atic errors (EMi); 2) random errors (EM 2 ), for which 
we assume that S c or 5' c can only be random errors of 
the Gaussian or the uniform type; 3) coexistence of both 
systematic and random errors (EM 3 ). In the next sec- 
tion, we shall present the results of numerical simulations 
and discuss the effects of imperfect gate operation on the 
DFT algorithm, and thus on the Shor's algorithm. 



III. INFLUENCE OF IMPERFECT GATE 
OPERATIONS 

We first discuss the influence of imperfect gate opera- 
tions in the initial preparation 

A^A^MO-O) = ^(|0) + |1) + 5x(|0) - |1») ® (|0) + |1) + 6 

If the errors are systematic, for instance, caused by the 
inaccurate calibration of the rotations, then 5\ = 62 = 
. . . = S n — 6. In this case, we can write the 2nd term as 
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where s stands for the number of l's, and 2s — n = s — 
(n — s) is the difference in the number of l's and 0's. 
Thus the results after the first procedure is 

2 L — 1 

lj2Qa)+5(2s-n)\a)) = ^=J2(l + 5 a )\a). (2) 

This implies that after the procedure, the amplitude of 
each state is no longer equal, but have slight difference. 
Combining the effect in the initialization and in the DFT, 
we have 

(1 + <5 a ) (1 + 5 c )e^ +5 > = (1 + 5")e^+ S >, 
where 5" = 5 C + 5 a . In the DFT, we have 
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where we have rewrite 5" as Sj here. Let P c denote the 
probability of getting the state |c) after we perform a 
measurement, we have 
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(1 + S' c )\~4,(^ E^ 1 + ^K 1 + cos[— r(m — k) + (mr + l)8' r< 
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From Eq. (3), we find that after the last measurement, 
each state can be extracted with a probability which is 
nonzero, and the offset I can't be eliminated. 

Eq. (3) is very complicated, so we will make some 
predigestions to discuss different error modes for conve- 
nience. Generally speaking, the influence of exponential 
error 5j is more remarkable than Sj, so we can omit the 
error Sj, thus 
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A. Case 1 

If only systematic errors (EMi) are considered, namely, 
all the (5j's are equal, then /(c) can be given analytically 
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The relative probability of finding c is 
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and if c = fc-, then 
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It can be easily seen that Hmj-^Q P c = -, which is just 
the case that no error is considered. 

When S takes certain values, say, 5 — |(fc — ^)vr where 
fc is an integer, then the summation in Eq. ([4J is on 
longer valid. In our simulation, 6 does not take these 
values. Here we consider the case where q — 2 7 = 128 
and r — 4. For comparisons, we have drawn the relative 
probability for obtaining state c in Fig.l. for this given 
example. We have found the following results: 
(i) When S is small, the errors do hardly influence the 
final result, for instance when c = fc-, then 
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the effect of different mode of random errors, we use two 
random number generators. One is the Gaussian mode 
and the other is the uniform mode. In this case, the er- 
ror has the form S = Sq + s, where So is the systematic 
error, s has a probability distribution with respect to 
c, depending on the uniform or the Gaussian distribu- 
tion. When 5q — 0, we have only random errors which is 
our error mode 2. When So =^ 0, we have error mode 3. 
For the uniform distribution, s ~ ±s ma x x u(0, 1) where 
m(0, 1) is evenly distributed in [0,1]. s ma x indicates the 
maximum deviation from So- For Gaussian distribution, 
s ~ N(Q, (To). Through the figure, we see the following: 

(1) When only random errors are present (So = 0), the 
peak positions are not affected by these random errors. 
However, different random error modes cause similar re- 
sults. The results for uniform random error mode are 
shown in Fig. 3. For the uniform distribution error mode, 
with increasing S max , the final probability distribution of 
the final results become irregular. In particular, when 
Smax is very large, all the patterns are destroyed and is 
hardly recognizable. Many unexpected small peaks ap- 
pear. For the Gaussian distribution error mode, as shown 
in Fig. 4, the influence of the error is more serious. This 
is because in Gaussian distribution, there is no cut-off 
of errors. Large errors can occur although their proba- 
bility is small. The influence of ao on the final results 
is also sensitive, because it determines the shape of the 
distribution. When ao increases, the final probability dis- 
tribution becomes very messy. A small change in ao can 
cause a big change in the final results. 

(2) When So 7^ 0, which corresponds to error mode 3, 
the effect is seen as to shift the positions of the peaks in 
addition to the influences of the random errors. 
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The probability distribution is almost identical to those 
without errors. 

(ii) Let us increase S gradually, from Fig. 2, we see that 
a gradual change in the probability distribution takes 
place. (Here, we again consider the relative probabilities) 
When S is increased to certain values, the positions of 
peaks change greatly. For instance at S — 0.05, there 
appears a peak at c=127, whereas it is P c = when no 
systematic errors are present. In general, the influence of 
systematic errors on the algorithm is a shift of the peak 
positions. This influences the final results directly. 



B. Case 2 

When both random errors and systematic errors are 
present, we add random errors to the simulation. To see 



IV. SUMMARY 



To summarize, we have analyzed the errors in Shor's 
factorization algorithm. It has been seen that the effect 
of the systematic errors is to shift the positions of the 
peaks, whereas the random errors change the shape of the 
probability distribution. For systematic errors, the shape 
of the distribution of the final results is hardly destroyed, 
though displaced. We can still use the result with several 
trial guesses to obtain the right results because the peak 
positions are shifted only slightly. However, the random 
errors are detrimental to the algorithm and should be 
reduced as much as possible. It is different from the 
case with Grover's algorithm where systematic errors are 
disastrous while random errors are less harmful pfl| . 
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FIG. 1: Relative probability for finding state c in the absence 
of errors. 
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FIG. 2: The same as Fig.l. with systematic errors. In sub- 
figures (1), (2), (3), (4), 8 are 0.02, 0.03, 0.05 respectively. In 
sub-figure (4), the curve with solid circles(with higher peaks) 
is the result with S = 0.1, and the one without solid cir- 
cles(with lower peaks) denotes the result with S = 0.33. 
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FIG. 3: The same as Fig.l. with uniform random errors. In 
sub-figures (1), (2), (3), (4), s max are set to 0.01, 0.03, 0.05, 
0.1 respectively. 
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FIG. 4: The same as Fig.l. with Gaussian random errors 
and systematic errors. In sub-figures (1), (2), and (3) r are 
set to 0.01, 0.03 and 0.05 respectively, and So = 0(without 
systematic errors). In sub- figure (4), both systematic and 
random Gaussian errors exist, where So = 0.33, r = 0.02. 



